Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 25 Febraary 2013 (MN MbX style file v2.2) 



Elemental abundances of intermediate age open cluster NGC 3680 



o 



\0' 

6 



in 

00 

o 

rn 

O 



X 



A. W. Mitschang ^^ f, G. M. De Silva ^ D. B. Zucker ^'^'^ 

^ Macquarie University Research Centre in Astronomy, Astrophysics & Astrophotonics 
^ Department of Physics & Astronomy, Macquarie University, NSW 2109 Australia 
^ Australian Astronomical Observatory, PO Box 296, Epping NSW 1710 



25 February 2013 



ABSTRACT 

We present a new abundance analysis of the intermediate age Galactic open cluster NGC 
3680, based on high resolution, high signal-to-noise VLT/UVES spectroscopic data. Several 
element abundances are presented for this cluster for the first time, but most notably we derive 
abundances for the light and heavy s-process elements Y, Ba, La, and Nd. The serendipitous 
measurement of the rare-earth r-process element Gd is also reported. This cluster exhibits a 
significant enhancement of Na in giants as compared to dwarfs, which may be a proxy for 
an O to Na anti-correlation as observed in Galactic globular clusters but not open clusters. 
We also observe a step-like enhancement of heavy s-process elements towards higher atomic 
number, contrary to expectations from AGB nucleosynthesis models, suggesting that the r- 
process played a significant role in the generation of both La and Nd in this cluster 
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1 INTRODUCTION 

Galactic open clusters (DCs) have long been established as im- 
portant fossils of both the dynamical and chemical evolution 
of the Galaxy. They are numerous, enable reasonably accurate 
dating, and eve n though th e ir life times are typically thought to 
be very short ( IJanes et al.l Il988h . significant numbers of old 
PCs have been detected ranging back to t he primordial Galaxy 
(|PheIpsetalill994l;lLiu& Chaboyer, 200^; iRandich et all |2009| ; 
iBragaglia et all l2008l : [Sestito et al.L l2008h . Hence they provide a 
time-resolved sample, across a range of Galacto-centri c radii, with 
which to probe the chemical evolution of the Galaxy Jpriel et all 
l2002h . Typically, [Fe/H] abundances are used as tracers of chem- 
ical evolution across OCs, and there are a large number of clus- 
ters for which these values have been computed, but only a mod- 
est sample of studies g o beyond this basic indicator (e.g., see 
ICarrera & PancinoluOl lb . into the realm of the a, r-, and s-process 
elements. 

These elements are now recognized to be important 
generational indicators, and form a basis by which to 
piece together the primordi a l buil ding blocks of the Galaxy 
JFreeman & Bland-Hawthoml 120020 . It has been demonstrated 
that their patterns, which become imprinted in the photospheres of 
cluster members when they are bom, can be used to link dissolved 
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associations to their natal environments JDe Silva et al.L l200q . 

2007; Bubar& King, 2010). 

The nucleosynthesis of these elements is believed to hap- 
pen i n the very late phases of stellar evolution (IWallerstein et all 
Il997h : specifically, Asymptotic Giant Branch (AGB) phase stars 
are thought to be responsible for the enrichment of some s-process 
elements including Y, Ba, Nd, and La, while high mass stars {M > 
8Mq) that are fated for death via core-collapse supernova (SN II) 
are responsible for a range of the a-elements including Mg, Si, Ca, 
and Ti. Type la supemovae (SN la) generate large amounts of the 
Iron-peak elements including Fe, Ni, Cr, Co and Mn. Several au- 
thors have established trends for these elements amongst thin and 
thick Galactic disk field stars for a range of me tallicities and ages 
teensbv et all 12001 I2OO3I : iReddv et all l2003h . enabling them to 
make conjectures about the relative rates and timescales at which 
these types of events occurred. 

NGC 3680, an intermediate age (~1.5 Gyr) OC situated 
about 8 kpc from the Galactic center in the thin disk, has been 
extensivel y studied photometrically and to some extent spectro- 
scopicallv. lNordstroem et al.l (119970 established membership prob- 
abilities based on radial velocities (RVs) and proper motions, 
and though there has been some spread in the literature over 
cluster m ean [Fe/H] values, from ~— 0. 17 dex to ~0.09 dex 
fNordstroe iii et al.lll997UBruntt et al.l 1 19991: IPasguini et all 1200 ll : 



20091) . 



lAnthonv-Twarog et al.L I2009I) . the general consensus from recent 
studies is that this cluster is slightly metal poor at ~— 0.08 dex. 
Most spectroscopic studies have been targeted at understand- 
ing the nature of the Li-dip observed in this and other clusters 
JAnthonv-Twarog et al.Ll2009l ; |Pasquini et al.Ll200lh , and thus there 



© 0000 RAS 



Mitschang & De Silva & Zucker 



Table 1. Target data 



ARES vs. SPLOT EW 



No. 



RA 

(h:m:s) 



Dec 

(d:m:s) 



Teff 

(K) 



log? 
(cm s~^) 



(km/s) 



11 


11:25:29.2 


-43:15:48.0 


10.88 


5100 


13 


11:25:16.1 


-43:14:24.3 


10.78 


4950 


20 


11:25:26.2 


-43:11:24.2 


10.10 


5200 


26 


11:25:38.0 


-43:16:06.4 


10.92 


5100 


27 


11:25:41.9 


-43:17:07.0 


10.73 


5050 


34 


11:25:38.6 


-43:13:58.9 


10.60 


5100 


41 


11:25:48.4 


-43:09:52.7 


10.88 


5100 


44 


11:25:49.8 


-43:12:16.0 


9.98 


4800 



3.1 


1.60 


3.0 


1.69 


3.2 


1.61 


3.3 


1.73 


3.1 


1.78 


2.9 


0.94 


3.2 


1.76 


2.6 


2.10 



24 


11:25:34.3 


-43:15:21.6 


13.77 


6650 


4.8 


2.13 


35 


11:25:38.1 


-43:13:27.3 


13.07 


6800 


4.4 


2.20 


5 


11:25:18.4 


-43:16:25.0 


12.83 


7100 


4.1 


2.20 



is a lack of publishe d abun dances beyond Li and Fe, though 
lAnthonv-Twarog et al.l ( 120090 derived Si , Ca, a nd Ni abundances 
for their program stars, and Pace et al.l (120081) reported several 
metal abundances for two dwarf stars. 

In this paper we present an abundance analysis of 8 giants and 
3 dwarfs based on high resolution spectroscopic data, covering 14 
elements, including several s-process and one rare-earth r-process 
element. Gadolinium. In Section |2] we describe our program tar- 
gets and observations and describe the methods used in determining 
differential abundances and uncertainties. In Section |3] we discuss 
membership, general characteristics and abundance trends of NGC 
3680, and finally in Section|4]we summarise our findings. 



2 OBSERVATIONS & ANALYSIS 

The targets in t h is stu dy were selected from the catalog of 
iNordstroem et al.l l ll997l) , as members of NGC 3680, using the 
membership determination criteria therein. Observations were car- 
ried out using the FLAMES multi-object fibre system feeding the 
UVES spectrograph (R~45,000) at the VLT as a part of service 
program 072.B-0331B over 9 nights in 2004. Table [T] gives de- 
tails of the targets in the study, including stellar parameters as 
determined below, where the top part shows giants and the bot- 
tom part dwarfs. Individual objects were observed at several dis- 
tinct times; all spectroscopic analyses herein were performed on 
exposure-weighted averages of all available spectra for each ob- 
ject. Care was taken to ensure matching of the wavelength scale for 
the different epochs observed before combining. Typical signal-to- 
noise (S/N) was ~130 for the giants and ~90 for the dwarfs (star 
Nos. 5, 24, 35). 



2.1 Equivalent width measurements 

Absorption lin e equivalent widths (EWs) were measured using the 
ARE:^ code l lSousa et al.L 120071) for automatic equivalent width 
measurements. Briefly, the code works by first determining the con- 
tinuum level using an iterative point rejection technique, normalis- 
ing the data, locating separate absorption features by looking for 
crossings in the third differential of the data and then fitting Gaus- 
sian profiles to these features. There are several input parameters, 
the optimal values of which were determined by comparing hand 



http://www.astro.up.pt/~sousasag/ares/ 




Figure 1. Comparison between ARES and hand measured EWs for a ran- 
dom sample of lines. The solid line indicates a one to one relationship. 



measurements of EWs in a typical spectrum using the splot func- 
tion in IRAIo until the best agreement was achieved. Figure[T]illus- 
trates the performance of the ARES code on these high resolution 
data. The results of EW measurement were further checked by re- 
viewing plots of the continuum and Gaussian profile fits for each 
line in each target. Where there were anomalous measurements 
(e.g. continuum contaminated by a spurious cosmic ray) the ARES 
measurements were replaced by hand measurements of those fea- 
tures. Table|2|lists the atomic parameters A, EP, and log gf for each 
species along with measured EWs for each star. 



2.2 Stellar parameter determinations 

The stellar parameters log g, Tcff, and microturbulent velocity ^, 
form the basis required to translate EWs to abundances or to fit syn- 
thetic spectra, based on assumed abundances, to observed spectra. 
We determined th ese parameters using spec troscopic methods and a 
grid of ATLAS9 JCastelli & KuruczlHooi) model atmospheres in- 
terpolated to give a resolution of 50K in Tes and 0. 1 in log g. Using 
the MOOCfl l lSnedenl 1 19731) driver abfind to force-fit abundances 
to our measured EWs, stellar parameters were determined by hand 
for each star in the typical way: Toh by balancing Fe I abundances 
against excitation potential, ^ by balancing Fe I abundances against 
reduced EW, and finally log g by requiring ionization balance be- 
tween Fe I and Fe II. On average across our eleven targets 120 Fe I 
lines and 11 Fe II lines were used in the determination of stellar 
parameters. 



2.3 Abundance Analysis 

The atomic line data for Fe and the other EW species are a sub- 
set o f those used by Bensbv et alj l l2003r) , Ba data were taken 
from iMcWilliamI ( 1 19981) while La, Nd and Gd were obtained via 



^ http://iraf.noao.edu/ 

^ http ://w w w. as .utexas . edurchris/moog. html 
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Table 2. Equivalent Widths for program stai's 



A 


ion 


EP 


logg/ 


11 


13 


20 


26 


27 


34 


41 


44 


24 


35 


5 


A 


















mA 












4799.41 


Pel 


3.64 


-2.23 


68.25 


75.40 


67.43 


74.30 


73.38 


49.78 


74.28 


85.48 


19.58 


12.63 


7.30 


4802.88 


Pel 


3.64 


-1.51 


88.31 


92.50 


88.08 


91.32 


92.29 


64.86 


93.29 


94.64 


47.72 


40.99 


28.40 


4808.15 


Pel 


3.25 


-2.79 


62.40 


70.62 


58.60 


66.12 


68.93 


49.96 


62.75 


73.66 


15.20 




4.51 


4809.94 


Pel 


3.57 


-2.72 


45.76 


52.83 


45.26 


53.19 


50.75 


39.79 


48.97 


57.12 


10.84 






4835.87 


Pel 


4.10 


-1.50 


102.19 


113.20 


100.48 


113.92 


111.55 


83.96 


59.23 


102.64 


36.46 


25.23 


19.81 




5087 5088 

Wavelength (A) 



5089 



Figure 2. ARES automatic line fitting for a small region of one of these 
high-resolution spectra. Only lines used in our analysis are shown here, but 
the code fits all features in the region for accuracy, including the features to 
the left and right of the line group shown here. 



the VALD databas^ll JKupkaet alll 19991) . Abundances for the ele- 
ments, Fe, Zn, Ca, Na, Mg, Cr, Ni, Ti, Si, and Al were then force-fit 
using abfind to all measured EWs for a given species. Values rela- 
tive to solar (i.e., [X/Fe], solar values derived as described below) 
are shown in Table |4] along with separate cluster means for giants 
(top) and dwarfs (bottom). 

For the species Ba, La, Nd, Gd, whose lines were too blended 
or for which the effects of hyperfine splitting preclude simple EW 
measurements, synthetic spectra were generated and compared to 
observations using the MOOG driver syiith. A range of abundance 
values were selected and compared by eye until a suitable fit was 
identified. Spectrum synthesis was not performed on the three 
dwarfs because of the lower S/N. The results of synthesis are re- 
ported together with other abundances in Table [4] 

Recent work by ISantos et al.l ( 120090 derived spectral param- 
eters for 3 giants in common with those in this study using VLT 
UVES data at R~50,000. The Tog, log g, and ^ values ob- 
tained from our analysis are systematically higher by, on aver- 
age, 250K, 0.4, and 0.3 km/s, respectively. To address any pos- 
sible systematics in terms of atomic parameters, we repeated our 
stellar parame ter determination using Fe atomic data fr om the 
same fine fist JSousaet alll2008l) as in lSantos et"aLl ( |2009|) . Given 
the differing spectral range in our data, we were only able to 

* http://vald.astro.univie.ac.at/~vald/ 
© 0000 RAS, MNRAS 000, 000-000 



utihze 181/263 Fel lines and 24/36 Fell fines. With this anal- 
ysis we came out with parameters compatible with our previ- 
ous measurements. To f urther test ou r Fe a tomic parameters, in 
the same manner as in ISantos et al.l ( 120091) , we derived stellar 
parameters for the HARPS solar spectrum "Ganymede" again 
with the same procedure followed above, obtaining very close to 
expected values (Tcft=5850±50, log g=4.4±0.1, ^=0.91 ±0.1, 
[Fe/H]=0.02±0.08). Therefore, the differences in stellar parame- 
ters are likely due to differences in data and code used to derive the 
them. 

Systematic errors affecting one's ability to compare abun- 
dance results from different studies come from a variety of sources, 
including choices of model atmosphere grids, atomic data and more 
subtle effects such as automated algorithms for continuum determi- 
nation. In order to minimise the contribution of such systematics, 
all abundance values presented here are reported relative compared 
to solar measurements as described above, using the "Ganymede" 
spectrum. We derive solar values via EWs for Fe, Na, Mg, Al, Si, 
Ca, Ti, Cr, Ni, Zn, and Y of 7.54, 6.43, 7.79, 6.52, 7.53, 6.42, 4.92, 
5.73, 6.28, 4.69, and 2.17, respectively. For the synthesised ele- 
ments Ba, La, and Nd, our computed values are 2.23, 1.12, and 
1.45. 



2.4 Errors and Uncertainties 

The sources of uncertainty in abundance values are related to the 
measurement of EWs (including the automatic continuum determi- 
nation), the atomic parameters (including excitation potential and 
log gf), and the determination of stellar parameters (which in turn 
are derived from Fe EW measurements, atomic parameters and 
choice of atmospheric model). 

Uncertainties for EWs were estimated via a Monte-Carlo 
method, by inducin g variations in the m ost sensitive ARES parame- 
ter, rejt (as noted bv lSousa et alj200'fl) , between reasonable values, 
(0.990,0.999) for giants and (0.885,0.995) for dwarfs (the differ- 
ence in mean rejt values between dwarfs and giants being a func- 
tion of S/N). We collated measurement differences from the mean 
across all lines in a given species for all stars. The resulting distri- 
bution was approximately Gaussian with the zero point being the 
EW used in analysis; the 90% limit of this distribution (i.e., 90% of 
all differences lie within this value) was then used in determining 
the abundance sensitivity for each species for a representative giant 
and dwarf. 

Uncertainties were conservatively assumed to be at the mea- 
surement resolution for Tog and log g, 50K and 0.1 respectively, 
and the uncertainty in ^ was taken to be 0.1 km/s. Abundance sen- 
sitivities are reported in Table |3]for all identified sources of error, 
where 5L and 5T are the uncertainty on log g and Tcff respectively. 
The quadrature sum of all these sources is reported as the total un- 
certainty. Due to the marked difference in quality of spectra for 
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Table 3. Abundance uncertainty sensitivities 





Fe 


Na 


Mg 


Al 


Si 


Ca 


Ti 


Cr 


Ni 


Zn 


Y 


Ba 


La 


Nd 


Gd 


5L±0.1 
5T±50 
5€±0.1 
5EW 


±0.01 
±0.04 
=F0.04 
±0.06 


±0.05 
TO.04 
=F0.03 
±0.09 


±0.01 
±0.03 
TO.03 
±0.06 


±0.01 
±0.03 
±0.01 
±0.06 


±0.02 
TO.01 

=fo.oi 

±0.07 


TO.01 
±0.04 
TO.04 
±0.07 


±0.02 
±0.04 
±0.03 
±0.03 


±0.01 
±0.04 
TO.02 
±0.05 


±0.01 
±0.03 
TO.04 
±0.06 


±0.03 
TO.01 
=F0.06 
±0.00 


TO.01 
±0.06 
TO.07 
±0.03 


±0.05 
±0.01 
=F0.03 
±0.04 


±0.04 
±0.01 
=F0.02 
±0.04 


±0.04 
±0.01 
TO.04 
±0.06 


±0.04 
TO.01 
=F0.03 
±0.06 


Totgiant 


±0.08 


±0.11 


±0.07 


±0.07 


±0.08 


±0.09 


±0.06 


±0.07 


±0.08 


±0.06 


±0.10 


±0.07 


±0.06 


±0.08 


±0.08 


5L±0.1 
5T±50 
5^±0.1 
5EW 


±0.00 
±0.03 

to.oi 

±0.07 


±0.01 
±0.02 
±0.01 
±0.12 


TO.01 
±0.02 
^0.01 
±0.04 




±0.00 
±0.01 

=fo.oi 

±0.07 


TO.01 
±0.03 
TO.01 
±0.06 


±0.00 
±0.04 

=fo.oi 

±0.06 


±0.00 
±0.03 
±0.00 
±0.16 


±0.00 
±0.03 
TO.01 
±0.07 


±0.01 
±0.03 

=fo.oi 

±0.00 


±0.01 
±0.04 

to.oi 

±0.01 










Totdwarf 


±0.07 


±0.12 


±0.05 




±0.08 


±0.06 


±0.07 


±0.17 


±0.07 


±0.03 


±0.04 











Table 4. Differential [X/Fe] abundances 



No. [Fe/H] 



Na 



Me 



Al 



Si 



Ca 



Ti 



Cr 



Ni 



Zn 



Ba 



La 



Nd 



Gd 



11 


0.04 


0.05 


-0.05 


-0.06 


0.04 


-0.15 


0.18 


0.01 


-0.07 


-0.34 


0.25 


0.16 


0.36 


0.56 


-0.39 


13 


-0.01 


0.13 


0.00 


0.12 


0.07 


-0.09 


0.22 


0.09 


0.01 


-0.38 


0.23 


0.16 


0.31 


0.56 


-0.34 


20 


0.06 


0.22 


-0.14 


0.08 


0.10 


-0.09 


0.24 


0.04 


-0.05 


-0.31 


0.10 


0.18 


0.24 


0.44 


-0.46 


26 


0.09 


0.12 


-0.01 


0.16 


0.01 


-0.03 


0.29 


0.11 


-0.03 


-0.43 


0.13 


0.12 


0.29 


0.51 


-0.49 


27 


0.04 


0.16 


-0.01 


0.10 


-0.02 


-0.11 


0.21 


0.10 


-0.10 


-0.50 


0.09 


0.06 


0.26 


0.56 


-0.34 


34 


-0.11 


0.04 


-0.15 


0.02 


-0.03 


-0.09 


0.07 


0.07 


-0.04 


-0.48 


0.19 


0.21 




0.21 


-0.64 


41 


0.03 


0.11 


-0.01 


0.18 


0.03 


-0.06 


0.25 


0.12 


0.00 


-0.42 


0.11 


0.13 


0.27 


0.55 


-0.43 


44 


-0.13 


0.18 


0.05 


0.27 


0.15 


-0.03 


0.21 


0.14 


-0.03 


-0.49 


0.24 


0.03 


0.43 


1.09 


-0.49 



Mean 


0.00 


0.13 


-0.04 


0.11 


0.04 


-0.08 


0.21 


0.09 


-0.04 


-0.42 


0.17 


0.13 


0.31 


0.56 


-0.45 


24 

35 

5 


-0.07 
-0.15 
-0.17 


-0.16 
-0.06 
-0.15 


-0.16 
-0.11 
-0.11 




0.00 
0.01 
0.01 


-0.05 
-0.04 
-0.01 


0.15 
0.07 
0.11 


0.04 
-0.07 
-0.16 


-0.03 
-0.10 
-0.09 


-0.29 
-0.23 
-0.16 


0.11 
0.12 
0.11 










Mean 


-0.16 


-0.10 


-0.11 




0.01 


-0.03 


0.09 


-0.12 


-0.10 


-0.20 


0.12 











our giants and dwarfs, we have computed the uncertainties for each 
separately. 



3 DISCUSSION 

3.1 Cluster membership and properties 

Although cluster membership has already been established for 
our target sample, it is prudent to confirm this, and in doing so 
also confirm our determination of stellar parameters. To this end, 
we performed a least squares fit to a set of Yale-Yonsei (Y^) 
JDemargue et al.l |2004|) isochrones against Tea and log g. Figure 
[5] shows the best fit isochrone, which yields a cluster age of 1.4 
Gyr, in good agreement with photometrically determined ages in 
previous studies of between 1.45 an d 1.75 Gyr (Nordstroem et al.l 
ll997l:lAnthonv-Twarog et al.ll2009h 

There is a clear outlier in Figure |3] the coolest of the dwarf 
stars (target No. 24), which has an anomalously high log g. The 
[Fe/H] for this star, though within the limits of the typical spread 
amongst all targets, is clearly deviant from the other dwarf stars. A 
review of the Fe lines used in determining the log g value revealed 
no unusual trends. The radial velocity for this target, at 2.44 km/s 
from the c luster mean, is on the out er envelope of those considered 
members dNordstroem et al.l Il997l) . raising the possibility that it 
has been erroneously classified as a member. We cannot dismiss 
the possibility that the isochrone placement of this target is due to 
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Figure 3. Best fit Y'^ isochrone to derived stellar parameters. The solid line 
represents the best fit isochrone of age 1.4 Gyr and metallicity —0.08 dex, 
while the dashed line shows solar level metallicity at the same age and dash- 
dot shows an isochrone of 1.3 Gyr at metallicity —0.08 dex for comparison. 
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NGC 3680: element abundances 



binaritv uNordstroem et alj ( 119971) marked it as a single line binary 
("SB 1" in their Table 1), however we find no such indications in our 
spectra, which consist of 8 epochs observed over 16 days. Given 
its position in the isochrone, along with the slight extremity of its 
radial velocity, we are led to reject it as a non-member based on our 
spectroscopic analysis. 

On examination of Table |4l we see two giants that strongly 
deviate in [Fe/H] from the rest of the giant population, specifically 
Nos. 34 and 44 at '-^0. 15 dex below the mean. Both of these stars 
also exhibit proportionally low Ba abundances, while star No. 44 
appears to have anomalously high Nd, and enhancement in several 
other elements (Mg, Al, Si, and La). This contrasts starkly with star 
No. 34's relative deficiency in Nd, Gd, Mg, Na and Ti. 

These two targets are squarely within the radial velocity cri- 
teria for membership and both have high membersh ip probabilities 
based on proper motion ( INordstroem et al.l 119971) . Both likewise 
do not appear anomalous on inspection of Figure|3] so it is unlikely 
the unusual abundance patterns are explained by mistaken mem- 
bership assignment. Binarity for No. 34 has bee n suggested based 
on a blue excess in photometric data tMenTiillio(L.et al., 1995), but 
it exhibits no indication of binarity in its spectrum. No. 44 is not a 
known binary system, ho wever we note that the bi nary fraction in 
this cluster is fairly high JNordstroem et al.l 1 1997b . Indeed, of the 
1 1 targets in this study five (Nos 11, 20, 27, 24, and 5) are consid- 
ered spectroscopic binaries, though they all clump quite tightly in 
Figure [3] 

3.2 Metallicity 

Ignoring the previously mentioned probable non-member, star No. 
24, and averaging between both dwarfs and giants as a single 
group, we obtain a cluster metallicity of [Fe/H]=— 0.03±0.02. If 
we treat dwarfs and giants as two separate but equally significant 
groups, we obtain the s li ghtly lower value of [Fe/H] =— 0.08±0.03. 
lAnthonv-Twarog et al.l | |2009|) . using a weighting scheme based 
on standard errors, report the cluster mean [Fe/H]=— 0.08±0.02, 
in perfect agreement with the second value above and only 
slightly poorer than the first. A perplexing result is that 
the relation between group means of giants and dwarfs, if 
considered sepa rately, is reversed in our analysis as com- 
pared to that in lAnthonv-Twarog et al.l 1 120091) . Where we com- 
pute [Fe/H]G=0.00±0.03, they obtain [Fe/H]G=-0.17±0.08, 
and while our dwarfs have [Fe/H]D=— 0. 16±0.05, theirs have 
[Fe/H]D=-0.04±0.11, differences of 0.16 and 0.12 dex respec- 
tively in opposite directions. It is important to note that these mean 
abundances for dwarfs are actually compatible within uncertain- 
ties, and that our value is based on only two dwarfs as compared 
to their thirteen. It is possible that the discrepancy between giants 
in our analyses can be attributed t o differences in stellar param- 
eter determination as described in lAnthonv-Twarog et al.l (l2009f) ; 
here we use ionisation balance, while they use photometric meth- 
ods. Referring to T able |3] had we adopted the stellar parameters of 
I Santos et alj ( l2009f ). the mean cluster metalicity would increase by 
~0.12dexto[Fe/H]=0.04. 

In spite of the differences in analysis mentioned above, and 
noting the paucity of abundances in the literature for this cluster, 
we now attempt to make comp arisons of other element abundances. 
lAnthonv-Twarog et al.l ( |2009|) were able to derive abundances of Si 
and Ni for both giants and dwarfs, and additionally Ca for their 
dwarfs. Between our studies for giants, Ni is compatible at close 
to solar values, yet Si in their analysis is significantly enhanced 
compared to ours, at [Si/Fe]G=0.22±0. 14, though this could be 



explained by the difference in numbers of lines used (their one 
line versus ten in the present analysis). For the dwarfs, both Si and 
Ni are at similar levels while Ca is 0. 1 dex below our m ean but 
still within uncertainties. Another study. IPace et al.l ( 120081) , derived 
abundances based on two dwarfs for Na, Al, Si, Ca, Ti, Cr and 
Ni obtaining 0.01±0.08, -0.08±0.04, -0.01±0.04, 0.04±0.06, 
0.04±0.08, 0.01 ±0.04, and -0.05±0.03 respectively. With the ex- 
ception of Na and Cr, which seem at odds, we see no major differ- 
ences amongst dwarf stars. Both Na and Cr for dwarfs exhibit the 
largest uncertainties amongst our entire stellar sample. 



3.3 Gadolinium 

While synthesizing Nd we noticed a feature in the spectra that could 
only be properly reproduced by modulating the abundance levels of 
the r-process element Gadolinium. This feature also appears and is 
well fit with our atomic data in the solar spectrum. To our knowl- 
edge, there exist no abundance measurement s of Gd in open clus- 
ters in the literature. ISneden & Parthasarathvl(ll983h . using several 
Gd features at 3549A and 3768A in the spectrum of HD 122563, 
obtained a simila r ly defi cient abundance of [Gd/Fe] =— 0.50, and 
lOen Hartog et al.l ( 120061) addressed the lack of measurements by 
improving abundance measurements of the Sun and deriving abun- 
dances for three metal poor giants obtaining values of —0.14, 
—0.42, and —1.08 respectively. Unfortunately, as these are not 
known cluster members and were selected for their low metal abun- 
dances, it is difficult to place those results in the context of the cur- 
rent study. Our measurements are based on a single line at 4463A, 
which is evident in the solar spectrum and likewise can be fit in the 
NGC 3680 giant spectra. 

Figure |4| shows synthesis results for Nd and Gd for both the 
Ganymede solar spectrum and a representative giant spectrum from 
our analysis; note the difference in line shape between the giant 
star and the Sun, which is likely caused by broader profiles of the 
Fe features on either side. Our computed solar abundance of Gd, 
using line data o btained from VALD, is 2 .77 which is well above 
that computed bv lOen Hartog et al.l ( |2006|) of 1.11. Effects such as 
hyper-fine structure or isotopic ratios may effect the computations 
significantly and were not included in our synthesis, as we are not 
aware of any analysis specifically involving the line at 4463 A. With 
the cluster giants at near solar metallicity, the robustness of our dif- 
ferential values seems enhanced; however effects due to tempera- 
ture and log g (i.e., comparing giants and dwarfs) may play a sig- 
nificant role in the shaping the Gd profile, making such comparison 
difficult. 

Eu synthesis proved to be too difficult for this cluster so Gd 
is the only pure r-process element we were able to measure. Gd 
may act as a proxy for the trends of other rare earth abundances 
in this cluster; however given the lack of specific knowledge of 
this species, in addition to our measured solar value well beyond 
literature solar values, our reported abundance must be taken with 
caution. 



3.4 Na enhancement in giants vs. dwarfs 

There is a clear divide between giants and dwarfs in terms of 
[Na/Fe] abundances, of order 0.25 dex, as seen in Table|4] Assum- 
ing that the stars in this clust er formed out of a well mixed com - 
mon natal environment (e.g. iFreeman & Bland-HawthomI 120020 . 
it is possible that the initial abundances of Na in RGB phase 
stars are modified and enhanced via the Ne-Na cycle, after which 
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Figure 4. Nd and Gd synthesis for the HARPS Ganymede and one of our 
target giant star spectra (plus and x symbols respectively). The thick red 
lines indicate our accepted fits, while the thin solid and thin dashed lines 
give an idea of the sensitivity, representing — 0.4/+0.5 dex for Nd and 
-1.3/+0.2 for Gd in star No. 13, and -1.6/+0.4 for Gd in the sun 



deep convect ion brings the newly formed Na to the photosphere 
JSalaris et al.l |2p02). It is unclear whether the temperatures re- 
quired for the Ne-Na cycle are reached in these stars, and thus a 
primordial explanation may be required. The commonly observed 
0-Na anti-correlation in Galactic Globular Clusters (GCs) may 
be related to this action, as the CNO cycle is thought to require 
the same temperatures/depths as the Ne-Na cycle, and is respon- 
sible for the depletion of O. Again, whether this observation is 
explained by stellar evolution or primordial abundance variations 
is still an open question. It is interesting to note, though, that the 
same anti-correlation has not been observed in OCs in the Galaxy 
JDe Silva et"aill2009l) . 

It is also possible that the abundance differences between 
dwarfs and giants are the result of non-LTE effects and no intrinsic 
differences exist. Na is particularly sensitive to non-LTE effects, 
though the m ost notable correct i ons ar e required at extremely low 
metallicities JAndrievskv et al.L |2007| : iBaumueller et al.L Il998h . 
Even after computing Na abundances using non-LTE models, 
lAndrievskv et al.l 1 120071) observe a Na enhancement in giants, also 
present in their LTE analysis, which they attribute to a mixing 
mechanism as above. 



3.5 The s-process element abundances 

Due to the lack of a range of s-process abundance measurements 
for NGC 3680 in the literature, we take this opportunity to discuss 
abundance trends for the elements Y, Ba, La, and Nd specifically. 
Figure|5]places the heavy element abundances of NGC 3680 within 
the context of other OCs and Galactic field stars w here possible. 
The th ick and thin disk samples were taken from iBensbv et al.l 
( 120051) while the cluster sampl e was compiled from references 
within lCarrera & Pancinol (1201 ID , using the mean abundance value 
from available high resolution studies for a given cluster NGC 
3680 exhibits higher than average La and Nd abundances, although 
several other clusters show more enhanced abundances for these 
two species. IC 4756 and NGC 2420 exhibit higher Nd abundances 
while NGC 2141 and Tombaugh 2 are more enhanced in La. NGC 
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Figure 5. [Fe/H] vs. heavy metal [X/Fe] abundances in a selection of open 
clusters and Galactic field stars. Filled squares represent cluster means as 
described in the text, while open and filled circles represent thin and thick 
disk field stars, respectively. NGC 3680 is shown by the encircled x sym- 
bol. Vertical and horizontal dashed lines mark solar abundance levels. The 
dotted line in each panel is a linear regression to the cluster abundance data, 
for illustrative purposes. 



2420 has both the lowest La and highest Y in the sample, and is 
additionally the most metal-poor cluster considered here. 

The two heavy s-process elements, La and Nd, appear to be 
more enhanced with respect to Fe as metallicity decreases amongst 
the cluster population. This trend is broadly consistent with that 
predicted in AGB stars due the the increase i n neutron exposur e 
with a decrease of the iron seed nuclei (e.g. see lBusso et al.ll200l[) . 
Ba does not seem to exhibi t such a trend, though the cluster abun- 
dance scatter is quite high. iD'Orazi et al.l ( 120091) also find no cor- 
relation in [Ba/Fe] with metallicity amongst clusters, but notice a 
significant anti-correlation with age, which they argue is evidence 
of a greater extra-mixing efficiency for lower mass stars producing 
the ^^C neutron source. NGC 3680, in our analysis, falls roughly 
0.06 dex b elow the group of in termediate age clusters shown in 
Figure 2 of lD'Orazi et al.l ( 120091) roughlv 0.06 dex below the group 
of intermediate ages, but still quite clearly within the uncertainties 
and thus does not challenge the age-Ba anti-correlation observed. 

The heavy s-process elements Ba, La, and Nd in the giant pop- 
ulation of NGC 3680 have increasing abundances with increasing 
atomic number, by roughly 0.2 dex at each step. Other OCs ex- 
hibit a similar monotonically increasing trend, as seen in Figure |6] 
which plots atomic number vs. mean cluster abundance for a sam- 
ple of clusters from Figure |5]which had measurements for all three 
elements, though none with a slope as consistently high as NGC 
3680. Since the s-process alone from models of AGB stars yields a 
clear peak at Ba res ulting in lower abu ndances for the other heavy 
s-process elements dBusso et al.ll 19991) across a range of metallic- 
ity and mass, another process is needed to explain the observed en- 
hancements. The most likely such process is the rapid neutron cap- 
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Figure 6. Heavy element trends with respect to atomic number for a subset 
of the cluster sample shown in Figure |5] NGC 3680 indicated by the red 
line and encircled x symbol. 



ture (r-) process, which quickly builds very large nuclei that decay 
over time to the stable isotopes beyond atomic numbers of around 
55. This process is thought to require high neutron densities, likely 
only available in the shocks of core-collapse supemovae. The Ba 
peak from AGB mod els is observed in so lar photospheric and me- 
teoritic data (e.g. see lAsplund et al.ll2009r) . implying a higher ratio 
of high-mass stars (those fated for the SN II) in the cluster's pro- 
genitor population, as compared to that of t he Sun. Other cluster s 
exemplify the Ba-peak (e.g. see Figure 2 of IPe Silva et alj|2009f) . 
indicating a range of formation scenarios that can be explored using 
abundance trends within the s-process elements. 

The sample in Figure |6] has a limited range in [Fe/H], owing 
in part to the limited availability of literature abundances for these 
heavy elements. Both the r- and s-processes have a dependence on 
metallicity due to the seeding nature of the Fe nuclei; more data are 
clearly needed at a range of metallicities lower than [Fe/H] ~— 0.5 
to explore this effect in more detail. 



4 CONCLUSIONS 

We have performed a detailed spectroscopic abundance analysis on 
the intermediate aged open cluster NGC 3680, obtaining differen- 
tial abundances of Na, Mg, Si, Ca, Ti, Cr, Ni, Zn, and Y for 8 giant 
and 3 dwarf stars, and Al, Ba, La, Nd, and Gd for the 8 giant stars. 
We find a combined metallicity of [Fe/H] =—0.03 dex, and for the 
giants in this study [Fe/H]=— 0.08 dex in good agreement with lit- 
erature values. Based on our analysis, we propose that the member- 
ship assignment of a single star (No. 24 in this study) be rejected 
both due to its position on the CMD and its anomalous metallicity. 

The serendipitous measurement of Gd in this cluster repre- 
sents the only measurement of a pure r-process element for this 
cluster, and likewise the only measurement of Gd for an open clus- 
ter that we are aware of. This may have important implications with 
respect to the progenitor population of NGC 3680, however due to 
measurement uncertainties it is difficult to comment on its value. 

A significant divide in [Na/Fe] abundances between dwarfs 
and giants is observed. This may be due to non-LTE effects which 
are difficult to quantify, but there is also the possibility that con- 
vection in the giants is bringing fresh Na to the photosphere, thus 
enhancing the [Na/Fe] abundance. If that is the case, this action 



may be related to the O-Na anti-coiTelation observed in Galactic 
GCs, but not OCs, though the mechanism responsible for this is 
still not well understood. 

An important result of this analysis is the measurement of a 
range of s-process elements (Y, Ba, La, Nd), for which measure- 
ments in open clusters remain limited in the literature. Compar- 
ing these abundances with available open cluster data, we find that 
NGC 3680 fits within the typical spread, with La and Nd having 
above average values. The Ba abundance from our analysis is con- 
sistent within uncertainties, though some what low for the cluste r's 
age, in the context of the recent finding bv lD'Orazi et al.l ( 120091) of 
an anti-correlation of [Ba/Fe] with age. 

We observe a step-like enhancement within the heavy s- 
process peak elements in this and other clusters. We argue that this 
observed trend indicates a substantial amount of r-processing con- 
tributed to the generation of La and Nd, likely in core-collapse su- 
pemovae in the progenitor populations of these clusters. 
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